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12 


energy  ci  teeing  fluids,  such  as  Minds  and  ccean  currents; 
(2)  the  potential  energy  of  tides  and  rivers;  (3)  the  heat 
generated  inside  the  earth;  (4)  the  direct  conversion  of 
sclar  energy  into  electricity  and  heat;  and  (5)  the  use  of 
sun  Maraed  Mater  in  conjunction  with  a source  of  cclder 
Mater  tc  provide  the  teaperature  differential  to  run  a vapor 
pewer  cycle.  Ncne  cf  these  sources  are  truly  Hfree".  lihat 
is  happening  is  the  trading  off  cf  the  transportation, 
processing,  and  ec vironaental  costs  of  a concentrated  energy 
source  fer  the  capital,  operating,  and  social  costs  of 
converting  a very  diffuse  source  into  a acre  concentrated 
fera  such  as  electricity. 

Xhe  earth  can  be  thought  of  as  a giant  heat  engine, 
abscrfcing  energy  frea  the  sun  and  reflecting  energy  as 
thcraal  radiaticn  tack  to  space  as  pictured  in  Bef.  2. 
Since  the  pcles  receive  less  energy  than  the  eguator,  the 
ataespbere  and  ocean  atteapt  tc  distribute  the  energy  aore 
evenly  ever  the  earth.  The  air  is  heated  ty  absorption  of 
sunlight  and  ty  contact  vith  the  surface  of  the  earth.  The 
heated  air  rises  thus  setting  up  surface  Mind  currents  as 
cccler  air  tries  tc  replace  the  rising  air.  The  shear 
stress  created  by  the  relative  aoticn  between  the  air  and 
the  surface  Mater  causes  the  Mater  to  aove  thus  creating  the 
surface  Mater  currents.  flany  currents  in  the  oceans, 
regardless  cf  depth,  originate  with  the  wind  shear  stress  at 
the  air-surface  interface.  The  circulaticn  patterns  cf  the 
oceans  beccae  extremely  complex  due  to  the  influence  of  the 
rotation  cf  the  earth  and  the  shapes  of  the  ccean  basins. 

In  general  the  surface  eaters  near  the  equator  are 
waned  by  the  sun  and  flew  towards  the  poles  giving  up 
energy  by  evaperation,  and  by  convccticn  tc  the  air  ard  by 
radiaticn  tack  to  space.  is  these  currents  cool,  scic  of 
the  water  becomes  dense  enough  tc  sink  and  then  flows  back 
to  the  eguater  alcng  the  bottoa. 
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The  relationship  between  salinity  and  temperature  and 
the  density  cf  seawater  is  such  as  tc  create  a stable 
stratification  of  the  ocean  in  aany  areas.  A salinity 
decrease  cr  a temperature  increase  reduces  the  density.  The 
total  effect  is  to  create  a warm  surface  current  overlaying 
the  cold  water  returning  from  the  poles.  In  the  tropical 
waters  the  surface  temperature  nay  exceed  80°F  (2€.7°C) 
while  2,00C  to  6,000  ft.  (914  to  1828  m)  dewe  the 
temperature  lay  be  about  40°F  <4.4°C).  This  temperature 
difference  can  tc  used  to  rur  a man-made  heat  engine,  like 
the  Ocean  Thermal  Energy  Converion  (OTEC)  power  plant 
ccccept. 

t 

The  proper  choice  of  a site  is  critical  to  the  cost  of 
such  a heat  engine.  Since  the  temperature  difference  is  low 
compared  to  conventional  thermal  power  generating  methods 
which  typically  use  temperature  differentials  exceeding 
50C°F  (26C°C ) , the  effect  of  a 1®F  (0.56°C)  loss  in  the  AT 
is  much  tore  costly  at  the  low  thermal  differentials,  which 
are  usually  less  than  50°F  (27.8°C)  in  the  ccean. 
Therefore,  the  vertical  temperature  profile  of  the  cceans 
will  ccntrcl  the  location  of  an  OTEC  power  plant.  The 
surface  temperatures  vary  with  location  and  with  time.  At 
most  places  in  the  ocean  there  is  a seasonal  change  of  the 
surface  water  temperature  of  5 to  8°F  (2.8  to  4.4°C)  except 
in  the  tropical  areas  that  remain  relatively  constant 
according  to  Bef.  3.  However,  all  regions  of  the  oceans  are 
affected  hy  wind-generated  mixing  of  the  surface  water  that 
changes  the  surface  temperature  by  a degree  or  so  and  a 
daily  ctcle  cf  a few  tenths  cf  a degree.  Some  cf  the 
possible  0.£.  near-shore  sites  being  considered  are  the 
waters  arcund  the  Eawaian  Islands  and  the  Gulf  Stream  off 
the  southeastern  U.S. 

In  the  past,  the  engineering  analysis  and  the  cost 
estiaaticn  have  keen  conducted  as  if  they  were  separate 
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functions.  Id  a systea  Kith  a saall  teapenature 
differential,  such  as  this  one,  the  engineering  design  and 
eccncaic  analysis  should  be  linked  together  as  a single 
procedure  because  of  the  capital  intensive  nature  cf  the 
ocean  theraal  power  plant.  For  exaaple,  the  heat  transfer 
egtaticn  (£  - (J A AT)  is  coupled  to  the  cost  equation  for 
heat  exchangers  (I  = K An)  through  the  surface  area  of  the 
tubes.  fecause  of  the  saall  AT  involved  the  area  atst  be 
very  large  tc  supply  the  saae  Q.  The  other  coaponents  are 
all  linked  tccether  by  the  engineering  and  cost  equations  so 
that  the  prcblea  cannot  logically  be  broken  into  an 
eguipaent  by  equipaent  optiaization.  The  present  analysis 
considers  a clcsed  vapor  power  systea  operating  cs  the 
theraal  pctertial  available  froa  the  ocean.  The  eccncaic 
optiaizaticn  cf  the  systea  is  carried  out  to  shew  what 
inferaatien  can  be  gotten,  what  conclusions  can  be  reached 
fzca  sucb  analysis  and  where  research  effort  sbculd  be 
expended  tc  iapxcve  the  design  of  the  systea.  This  study 
does  not  atteipt  tc  predict  the  cost  of  sucb  a systea  cr  to 
aake  specific  recoiaendations  about  the  paraaeters. 

A scheiatic  diagraa  of  the  systea  is  shewn  in  Figure  1. 
The  wan  surface  water,  the  Gulf  Streaa  for  exaaple,  is 
potFed  threueb  a heat  exchanger.  There  scie  energy  is 
transferred  by  boiling  the  working  fluid  such  as  aaacnia. 
The  working  fluid  is  piped  to  a turbine  where  soae  of  the 
energy  is  converted  to  aechanical  energy  then  to  electrical 
energy  is  a generator.  The  unavailable  portion  of  the 
absorbed  energy  is  rejected  to  the  environaent  in  the 
condenser  and  the  working  fluid  is  returned  to  the  liquid 
state.  Tbe  lew  teaperature  in  the  condenser  is  aaintained  by 
paaping  cc Id  ocean  water  up  froa  the  depths.  Finally  the 
acxking  fluid  is  puaped  back  to  the  boiler. 

Tbe  Ccear  Ibertal  Energy  Conversion  progran  has  received 
Bach  attention  in  the  last  several  years  froa  aany 
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researchers  funded  at  first  primarily  by  the  Rational 
Science  Fcurdation  (NSF) . Since  the  establishment  cf  the 
Energy  fieseatch  and  Development  Administration  (EEC A)  , new 
research  ccitracts  are  funded  mainly  by  this  agency.  The 
C.S.  Kavy  is  alsc  funding  some  research  in  the  OT EC  area. 
Bef.  4,  5 and  6 give  a gccd  summary  of  the  work  that  has 
been  done,  bbile  much  of  the  work  has  concentrated  on  the 
heat  exchangers,  the  other  components  and  the  systei  as  a 
vhcle  are  receiving  attention.  Some  researchers  are 
concentrating  on  the  environmental  and  legal  guesticns  of 
such  systems,  ether  teams  are  considering  hew  to  use  the 
energy  supplied  from  an  OTEC  power  plant. 

In  ary  design,  trade-offs  must  be  made.  Capital  cost 
can  be  traded  for  operational  costs.  Pumping  power  can  be 
traded  fer  heat  exchanger  size.  Trade-offs  can  involve 
items  that  have  a common  measure  such  as  dollars;  scce  ether 
trade-offs  are  very  hard  to  measure  in  common  terns  such  as 
sccial  and  ec viroEnental  costs.  The  trade-effs  may  involve 
crly  one  ccapcnent  such  as  a pump  or  a turbine,  but  most 
often  the  trade-effs  affect  more  than  cne  compcnent  or 
pcssihly  the  entire  system. 

The  trade-eff  process  should  continue  until  an  optimum 
design  is  ettained.  From  Bef.  7,  "The  aim  cf  optimization 
is  the  selection,  cut  of  the  multiplicity  of  potential 
solutions,  of  that  solution  which  is  the  best  with  respect 
tc  some  well-defined  criterion.”  In  order  for  the  plant 
design  tc  oe  optimum,  the  components  must  be  cptirized 
within  the  context  of  the  entire  systes.  Some  cf  the 
parameters  fex  a piece  of  cguipaent  may  have  no  direct 
effect  cn  the  optimum  design  of  another  portico  cf  the 
system,  but  seme  parameters  will  have  sigoificant  effects. 
Pox  example,  a heat  exchanger  is  smda  up  cf  a large  ruaber 
of  tubes  cf  a given  diameter  and  length  and  is  designed  to 
have  a certain  heat  transfer  rate  for  a given  flow  cate  of 
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water  pumped  through  the  tubes.  The  diameters,  length  and 
the  number  cf  tubes  may  be  optimized  for  lowest  cost  cf  the 
heat  exchanger.  This  approach  leaves  out  the  etfect  that 
the  tube  diameter  and  the  tube  length  have  cn  the  pressure 
drop  across  the  heat  exchanger.  Minimizing  the  cost  cf  the 
heat  exchanger  involves  ainimizing  the  surface  area.  This 
may  lead  tc  increasing  the  pressure  drop  across  the 
exchanger.  Ike  higher  the  pressure  drop  the  more  powerful 
and  the  xcre  expensive  the  pump  must  be  to  maintain  the  £1oh 
velocity.  Ibis  is  just  one  example  cf  the  hundreds  of 
trade-offs  possible. 

Unfortunately,  what  happens  in  many  cases  is  tha  *•  each 
group  cf  experts  nerking  on  a particular  component  optimize 
that  itex  aicund  some  given  conditions  that  may  have  teen 
set  ty  teams  meriting  cn  some  other  section  cf  the  system. 
Little  analysis  is  made  of  horn  the  entire  system  resperds  to 
the  coupling  censtraints.  Ill  that  the  puip  manufacturer 
Hants  is  that  the  buyer  tell  him  what  the  flow  rate,  head 
and  service  conditions  are  and  he  mill  gucte  a price. 
Eefcre  the  buyer  can  decide  on  the  flow  rate  and  bead  he 
should  first  know  how  the  ccst  cf  the  pump  varies  as  flov 
rate  and  head  vary  so  that  those  tno  parameters  can  be 
optimized  in  the  context  of  the  system. 

The  problem  then  becomes  cne  of  aeguiring  sufficient 
irfermatien  cn  hek  costs  vary  as  certain  parameters  change. 
Each  manufactures  is  able  to  make  estimates  cf  costs  fer  his 
pzcduct  but  he  is  reluctant  to  devulge  his  information  to 
buyers  cr  competitors.  This  makes  constructing  cost  curves 
a hazardous  (recess  at  best,  since  the  data  must  be  gathered 
by  ether  me t heds  that  may  cause  large  unexplained  variances. 

Sometimes  the  optimization  analysis  is  made  uith 
assumptions  that  make  the  answers  cf  gjestionable  value, 
for  example,  sometimes  the  heat  transfer  coefficient,  U,  is 
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assumed  fixed  in  the  equation  q * UA  AT.  Eut,  as  shewn  in 
Fet.  8,  for  toiling  the  heat  transfer  coefficient  is  very 
sensitive  to  the  temperature  difference  with  the  result  that 
q = U A ' AT  * • 3 5 . The  The  accuracy  of  the  engineering  equations 
should  te  kept  in  Bind.  Most  of  the  equations  are 
correlations  tc  fit  experimantal  data.  The  correlation  was 
developed  for  scie  particular  set  of  data  and  aay  disagree 
signif icartly  when  compared  to  data  taken  by  ether 
experimenters . The  linitations  of  the  various  relations 
Bust  be  kept  in  Bind  wh«»n  reviewing  the  results.  For  a 
systeB  tbe  size  of  the  OTEC  plant,  a pilct  plant  should  be 
constructed  in  order  to  prove  cut  the  answers  given  by  the 
calculations  before  expensive  mistakes  are  Bade. 

As  a design  progresses,  the  analysis  begins  with  a 
general  picture  containing  many  simplifying  assumptions  and 
proceeds  fereward  with  Bore  details  added  at  each  stage. 
This  report  is  the  second  one  written  at  tne  Naval 
Pcstraduate  fcbccl  on  the  subject  of  OTEC.  Commander  Furman 
Sheppard,  CEB,  in  his  analysis.  Ref.  9,  attempted  tc  show 
what  kinds  of  information  could  be  gained  froa  a ccibined 
thermal  ecocciic  model.  He  advocated  the  use  cf  zone 
analysis,  as  developed  by  Wismer  in  Ref.  10,  for  large 
complex  systets.  In  Hismer's  approach  to  2cre  analysis,  the 
systeB  is  broken  down  into  zones  containing  one  cr  mere 
components.  The  zenes  are  connected  by  linking  variables. 
Two  methods  cf  optimization  are  possible.  For  the  first 
method,  in  each  zone,  the  zone  parameters  are  optimized 
while  the  lirking  variables  are  held  fixed.  Then,  the  zone 
parameters  arc  held  fixed  while  the  linking  variables  are 
optimized.  Tbc  process  iterates  between  optimizing  the  zone 
parameters  and  the  linking  variables  until  convergence 
criteria  are  satisfied.  In  the  second  method  the  zones  are 
cut  off  trci  each  other  and  each  linking  variable  is 
separated  into  two  variables,  one  on  each  side  cf  the 
cutting  plate.  Tbe  optimization  proceeds  until  the  linking 


variables  ct  each  side  of  the  cuts  are  equal. 


Sheppard's  systen  did  not  produce  vapor  in  a boiler  tut, 
instead,  aaictainec  the  working  fluid  in  the  liguid  state  in 
tte  wan  seawater  heat  exchanger  and  vaporized  the  working 
fluid,  aoicnia  in  his  nodel,  by  expanding  it  through  a 
"black  tcx"  called  a vaporizer.  Hxs  optimization  analysis 
was  confined  to  the  zone  containing  the  warn  seawater 
exchanger,  the  feed  punp,  and  the  warn  water  circulating 
f unp. 


£ . CrJKUVIS 


The  present  analysis  inprcves  the  previous  nodel 
analysis  ty  replacing  the  liguid-to-liguid  heat  exchanger 
anc  the  vaporizer  with  a boiler.  The  boiler  nodel  is  tc  be 
cf  sufficient  detail  so  that  a realistic  desigr.  is 
represented.  The  ether  xajor  conpcnents  of  the  systen,  such 
as  the  turbirc,  the  condenser,  and  the  condenser  circulating 
puips,  are  included.  A nonlinear  pregranning  technique  is 
tsec  tc  peifcin  the  thernal-econonic  optinization  of  the 
ccnpleted  systen.  The  research  reported  here  is  intended  to 
show  the  usefulness  of  optinization  analysis  in  the  casign 
of  a lew  thernal  difference  systen  and  tc  denonstratc  some 
cf  the  picbleis  that  nust  be  considered  in  this  type  of 
analysis. 
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II.  ISI  BODE  I 


A.  GENESIL  CESC&IETION 


Vapor  f c we  r generating  cycles  usually  ccnsist  of  four 
processes.  Energy  is  absorbed  fro*  a high  temperature 
source  and  is  used  to  vaporize  a working  fluid  in  a toiler. 
Next , werk  is  extracted  by  expanding  the  vaper  in  an  engine, 
exhaustinc  tie  vaper  at  a lower  teaperature  and  pressure. 
The  working  fluid  is  then  condensed  back  tc  the  liquid  state 
by  rejectirg  energy  to  the  lewer  teaperature  energy  sink  in 
the  third  prccess.  Finally,  the  liquid  is  returned  to  the 
higher  presscre  of  the  boiler  by  a puap.  In  the  fcllcwing 
analysis,  tie  energy  source  is  a current  of  ocean  water 
teated  by  tie  sun  and  the  energy  sink  is  assuaed  to  be 
cclder  ccean  water  from  deep  currents.  The  werking  flcid  is 
assuaed  tc  be  propane  although  aany  other  fluids  are 
possible,  such  as  aaacnia  or  the  Freons. 

Six  eexpenents  are  aodelcd  in  the  following  analysis. 
The  hot  side  heat  exchangers,  the  boilers,  vaperize  the 
working  fluid  which  is  expanded  through  a turbine  exhausting 
to  a condenser  that  rejects  the  unused  beat  tc  the  cold 
sixk.  Three  puaps  are  included,  a feed  puwp  returns  the 
werking  fluid  to  the  boiler  and  two  puaps  circulate  the  wara 
amc  cold  waters  through  the  boiler  and  the  condenser, 
respectively.  The  arrangement  of  the  systea  is  shown  in  Fig. 
2.  The  nedes  are  nuabered  to  correspond  with  the  fluid 
state  points  shown  later,  in  Fig.  3. 
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:IGURE  2.  SYSTEM  ARRANGEMENT 


The  analytical  Model  consists  of  the  equations 
describing  tie  system  and  a set  cf  parameters  to  be 
optimized  sc  as  tc  minimize  the  capital  cost  of  the  plant 
subject  tc  certain  constraints.  The  objective  function,  the 
function  tc  te  minimized,  is  the  sun  of  the  capital  costs  of 
the  pumps,  the  boiler  and  the  condenser.  Although  the 
turbine  is  a significant  capital  cost  and  mould  affect  the 
design  cf  the  system,  it  is  not  included  in  this  «odel. 
There  are  tnc  classes  cf  constraints.  The  first  class  is  the 
explicit  constraints  that  specify  the  upper  and  lower  bounds 
• on  the  parameters.  The  other  class  of  constraints  is  the 
implicit  constraints.  These  are  the  engineering  eguations 
that  describe  the  operation  of  the  power  plant.  For 
example,  a constraint  might  say  that  the  temperature  inside 
the  boiler  must  te  less  than  the  temperature  of  the  warm 
seawater  because  the  equations  fail  if  the  requirement  is 
net  satisfied.  Another  example  would  be  a power  balance 
that  requires  that  the  sum  of  the  energy  flews  into  and  out 
of  the  system  te  zero.  Additionally,  they  proscribe  the 
region  that  the  optimization  rcutine  may  search  within  for 
the  minimum  cf  the  problem.  The  implicit  constraints  are  of 
two  types,  equality  and  inequality  constraints.  Inequality 
constraints  set  one-sided  bounds  on  the  feasible  region. 
Seme  of  the  inequality  constraints  are  restrictions  placed 
on  the  maximum  or  minimum  values  of  the  parameters  dee  to 
judgements  made  cn  the  basis  of  practical  engineering 
consideration.  The  equality  constraints  are  relationships 
that  must  be  strictly  satisfied  at  the  optininum.  The  system 
of  egoatiens  describing  the  model  fit  intc  the  following 
fen: 


, 

■ 

jia. 
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niriiize 


f (X) 


subject  to:  g (X)  £ 0,  i=1,2, ....  ,m 
i 

h (X)  = 0,  i = m+1 , . . . . ,m+n 
i 

where  X is  tke  vector  of  the  parameters  tc  be  optiaized, 
f ( X)  is  the  capital  cost  function,  (X)  are  the  inequality 
constraints  and  hj  (X)  are  the  equality  ccnstraints.  The 
parameters  are  the  dimensions  of  the  kcilers  and  the 
condensers,  the  boiling  and  condensing  pressures,  the  sea 
water  velocities,  and  the  mass  flew  rate  cf  the  werking 
fluid  tkrcucb  the  system. 

The  kasis  fer  the  economic  framework  must  be  set  cut  at 
the  start.  ill  measurable  costs  are  relevant  to  the  final 
design  cf  tke  power  plant,  but,  for  some  types  of  analysis 
and  design,  tany  ccsts  can  be  ignored.  The  life  cycle  cost 
of  the  Cl E C must  be  considered  before  the  decision  is  made 
tc  commit  the  U.S.  and  the  rest  of  the  world  tc  the 
ccr segue r.ces  cf  chcosing  this  energy  generation  scheme  as 
even  a partial  solution  to  the  energy  problem.  The  life 
cycle  cost  consists  of  the  capital  costs,  the  financing 
costs,  and  tke  operating  costs.  The  trade-eff  of  operating 
cost  and  capital  cost  is  not  considered  in  this  research, 
because  the  operating  cost  is,  for  the  most  part, 
independent  cf  the  overall  dimensional  characteristics  of 
the  components.  Operating  costs  are  dominated  by  the 
structure,  site  location,  and  maintenance  reguireaents. 
Such  things  as  the  kind  of  instrumentation,  controls, 
bearings,  acd  structural  materials  have  mere  effect  cn  the 
operational  ccsts  than  the  length  of  the  toiler  or  the 
diameter  cf  the  toiler  shell.  Therefore,  cnly  the  capital 
ccsts  of  tke  majer  components,  except  the  turbine,  are 
considered. 


24 


Tc  be  useful  for  model  analysis  cost  data  rust  be 
transformed  into  a cost  estimating  relationship  (CEB)  that 
shews  hew  cost  varies  with  some  variable  or  some  comb i ration 
or  varieties.  The  costs  must  be  demonstrated  to  be  highly 
correlated  with  the  variable  chosen.  Since  the  cost 
relation  is  cnly  a correlation,  the  bounds  on  tne  accuracy 
of  the  cguaticc  should  be  specified  by  those  who  develop  it. 
Since  the  relation  is  arrived  at  from  a finite  body  cf  data, 
the  CEB  is  directly  applicable  only  in  the  range  of  the  data 
collected.  Eitrapolation  cf  the  CEB  outside  the  ratge  is 
cargercus  because  there  is  nc  data  in  that  region  tc  support 
the  assumed  curve.  However,  the  value  of  CER's  is  the  use 
of  them  tc  enable  the  estimator  to  make  a better  analysis  of 
a new  design  than  he  could  if  the  information  is  not 
available.  The  relevant  range  of  the  CEB  must  be  kept  in 
■ind  wher  analyzing  the  solutions  from  the  standpoint  cf  the 
reliability  cf  the  answers.  Another  factor  tc  consider  when 
using  cost  estimating  relationships  is  whether  the  equipment 
the  data  was  taken  from  is  comparable  to  the  system  the 
analyst  is  considering.  If  the  design,  the  technology,  or 
the  applicaticn  of  the  proposed  equipment  is  significantly 
different  ftem  that  from  which  the  data  is  taken,  the  CEB 
may  give  unreliable  answers. 

This  analysis  begins  with  the  development  of  the 
engineering  eguatiens  that  constrain  the  system.  This  is 
followed  by  the  development  of  a cost  framework  based  on 
capital  costs  as  the  function  to  be  minimized.  The 
engineering  and  economic  relations  are  linked  together  in 
the  objective  function  and  the  constraints.  A description 
of  the  algorithm  that  performs  the  optimization  is  included. 


B.  CYCLE  A B J 1 ISIS 
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The  tuilcing  of  the  model  begins  with  the  analysis  of 
the  thencdyramic  cycle.  Figure  3 is  a temperature-entropy 
diagram  cf  tfe  wcrking  fluid  on  which  the  cycle  is  shown. 
These  numbered  nodes  represent  the  thermodynamic  state 
points  at  that  location  in  the  system  and  correspond  tc  tne 
numbers  ir  Fig.  2.  The  boiling  and  condensing  processes  are 
assumed  tc  be  reversible  processes.  In  the  pump  arc  the 
turbine  the  reversible  process  is  indicated  by  the  subscript 
"s"  on  the  state  point  number.  In  the  actual  system  there 
are  irreversible  losses  inside  each  component  that  are  left 
out  of  tte  analytical  model  or  lumped  in  with  ether  losses. 
The  pressure  drops  due  to  frictional  losses  in  the  toiler 
and  the  condenser  are  neglected. 

State  point  1 assumes  that  the  liguid  leaving  the  toiler 
is  saturated.  State  points  2s  and  2 represent  the 
conditions  leaving  the  feed  pump.  If  there  were  no  losses  in 
the  pump,  the  fluid  would  be  at  2s.  The  fluid  enters  the 
bciler  tc  te  heated,  first,  to  saturation  conditions  at 
state  point  3 then  boiled  at  essentially  constant  pressure. 
State  point  4 shows  the  vapor  leaving  the  bciler  as  a 
saturated  vapci.  lhe  lines  from  point  4 to  5s  and  5 
represent  tte  expansion  through  the  turbine.  The  path  tc 
state  pcirt  5 accounts  for  the  losses  due  tc  friction  and 
ether  irreversible  processes  in  the  turbine.  The  remaining 
heat  is  then  rejected  from  pcint  5 to  point  1,  again  assumed 
to  be  a constant  pressure  process. 

The  ttecretical  heat  transfer  rate  in  the  boiler  and  the 
ccndenser  is  found  by  subtracting  the  enthalpy  at  state 
pcints  2 and  1 from  the  enthalpy  of  state  points  4 and  5, 
respectively,  and  multiplying  by  the  total  mass  flow  rate  of 
the  working  fluid. 


Q.b  * - H*) 


(1) 
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ENTROPY 

FIGURE  3.  THERMODYNAMIC  CYCLE 


^re  j = Hu/f  (H 5 ~ ) 


(2) 


The  enthalpy  at  state  point  5 is  deterained  fcy  first  finding 
the  acisture  content  of  the  fluid  exiting  the  turbine. 


sf  st 


The  enthalpy  at  pcint  5s  is  deterained  froi 


HSs  8 H,  ♦ i H, 


This  aliens  f5  to  te  found  fren 


* - ' HS* 

'lur  ~ 


H H ' »5 


where  the  efficiency  of  the  turbine  is  assuaed  knewr.  Th< 
pener  out  of  the  turbine  is  non  known. 


$tur  ~ ~ Hj)  (6) 

for  the  feed  puap  the  change  in  enthalpy  is  accounted  for  by 
the  change  ir  pressure  assuming  negligible  rise  it  the 
teaperature.  The  enthalpy  of  point  2 is  determined  fret 


V1  (P»s  - Pi ) 


H2  - H, 


where  the  efficiency  is  assuaed  known.  The  power  reguired  to 
drive  the  puip  is  deterained  from 


*fp  = M^<H2 


- «l) 


lc  gain  a better  feeling  fcr  the  low  thermal  potential 
systeis  and  tc  take  comparisons  with  other  systems  several 
efficiencies  are  worthwhile  calculating.  ifce  first  is  the 
Carnot  efficiency  which  defines  the  maximum  efficiency 
possible  for  a system  operating  between  two  temperature 
liiits  ir  decrees  absolute,  °H  (°K). 


n . ^eold 

rear  » 1 - (9) 

^"hot 

The  next  efficiency  calculation  takes  into  account  the 
theoretical  tnermcdynamic  cycle  and  is  the  ratioof  the  net 
power  out  cf  the  cycle  to  the  beat  taken  up  by  the  cycle 


(10) 


Ibis  acccrnts  for  the  thermodynamic  inefficiencies  .inherent 
in  the  cycle.  Finally,  the  system  efficiency  takes  into 
acccunt  the  ether  internal  consumers  cf  power:  the 
circulaticn  pomps,  the  pipe  losses,  and  any  other  lesses, 
such  as  the  antifouling  equipment.  This  calculation 
accounts  fcr  all  of  the  mechanical  and  electrical 
ine f £ icierci es  of  the  system. 


NET  PONES  AT  THE  EUS  t AS 

Uflb 


01) 


For  example,  if  the  hot  water  is  at  75°F  (23.9°C)  and 
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ccld 

water 

is  assuaed 

efficiency 

is  6.58 

. Per  one  typ 

shewn 

in 

Table 

5,  where  th 

assumed  tc 

be  65* 

efficient,  th 

and 

the  s 

ystea  ef 

ficiency  is  2 

at  40°F  (4 . 4°C)  the  Carnot 
ical  set  of  assumptions,  as 
e turbine  and  the  puaps  are 
e cycle  efficiency  is  2.6* 
.48  for  one  feasible  design. 


C.  ENGINEERING  f R 1HEU0BK 


This  section  describes  the  design  o i t. 
make  up  the  system.  The  level  of  deta 
prevides  a physically  realistic  aodel  and 
lie  pipe  friction  lesses  on  the  cold  and  h 
well  as  tte  pipe  lesses  in  the  working 
neglected.  lhe  turbine  is  included  o 
transducer  and  no  atteapt  is  aade  to  oktai 
diiensicES,  cr  fluid  flow  characteristics, 
of  the  turbine  and  the  pumps  are  assuaed  a 
Ike  beat  absorbed  in  heating  the  wor 
saturated  state  is  a small  portion  of  the 
up  in  the  toiler;  therefore,  the  sensi 
fluid  is  neglected. 


1*  U&  Iciler 

The  tciler  consists  cf  a nuaber 
heat  exchangers  connected  in  parallel  as 
Figure  4 shows  the  internal  arrangeaent  of 
tute  bundle  is  horizontal  and  boiling  occ 
surface  cf  tie  tubes.  The  tube  bundle  is 
working  fluid.  The  working  fluid  is  aainta 
Just  cover  tte  top  rew  of  tubes.  The  warm 
through  tte  inside  of  the  tubes. 


he  components  that 
il  of  the  analysis 
clcses  the  loop. 


ct  water 

inlets  as 

fluid 

side  are 

nly  as 

an  energy 

n costs. 

ph  ysical 

The  eff 

ic iencies 

s known  constants, 
king  fluid  to  the 
tctal  energy  taken 
ble  heating  cf  the 


cf  tube-in-shell 
shewn  in  Fic.  2. 
cne  shell.  The 
urs  on  the  outside 
submerged  in  the 
ined  at  a level  to 
seawater  is  pumped 
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Ife  clearance  between  the  tubes  is  fixed  at  0.5 
inches  (1.21  cb)  to  allow  for  flaring  of  the  tube  ends  or 
ether  space  and  strength  reguireaents.  Sioilarly,  the  tube 
thickness  is  fixed  at  0.028  inches  (0.071  ca) , since  no 
strength  reguireaents  due  to  differences  between  internal 
and  external  pressures  on  the  tubes  or  nateiial  loss  cue  to 
cccrosion  cf  the  tube  surfaces  are  accounted  for. 

It  the  tep  of  the  tube  tank  the  leaving  vapor 
certains  a large  aacunt  of  liquid  droplets  which  *ust  be 
eliminated  either  in  the  boiler  or  by  an  external  separator. 
In  the  acdel  the  separation  takes  place  in  the  teller  ny 
allowing  sufficient  space  above  the  tube  bank.  Falen  and 
Snail,  fief.  11,  presented  a aethed  for  findieg  the  necessary 
separation  vcluae  that  considers  the  effects  or  surface 
tersicn,  density  of  the  vapor  and  the  liquid,  and  the  Bass 
flew  cf  tie  working  fluid.  First  the  allowable  vapor  lead  is 
found  fret 

V L * 2290A  F 2 ]°‘5  (12) 

^ 16.86x10-*  - p3  )J 

Dapcr  load  is  an  expression  for  the  aass  flew  rate  cf  vapor 
tnreugh  a urit  vcluae.  The  expression  for  the  naxiaua  vapor 
lead  Halts  the  flew  rate  so  that  the  aoisture  droplets  are 
able  to  settle  out  of  the  vapor  before  it  leaves  the  bciler. 
lkcc  the  required  seqaental  area  is  found  fret 


S A 


VL  I S 


(13) 


The  segaental  jc*a  is  the  cross  sectional  area  froa  the  top 
of  the  tube  bank  to  the  top  of  the  shell.  Kncwing  the 
segaental  area  the  height  of  the  tube  bank  above  the  Bain 


diaaeter  is  fcund  froa 


2 (2  hs\  2 2 0.5 

S A = C .25  Es  ABCOsf  — J-  hs(0.25  Es  - hs  ) < 1 4> 


Falsi)  a cd  Stall  also  present  a aucb  simpler  aethod  of 
deterainirg  tbe  total  height  of  the  tube  bank.  The  height 
of  the  separation  volume  is  set  at  40X  cf  shell  diaaeter 
which  results  in  "bsM  being  egual  to  0.1  Ds.  Either  eethod 
car  be  used  aloce  or  the  first  method  can  be  usee  as  a 
constraint  or  the  seccnd.  For  this  problem  the  segiental 
area  is  fcund  by  both  methods  with  the  first  used  as  a check 
cn  the  values  given  by  the  second. 

Sines  the  heat  flux  calculation  is  based  on  a single 
tube,  the  nuiber  of  tubes  per  shell  aust  be  determined  in 
order  tc  find  the  total  energy  absorbed  by  the  heat 
exchanger.  The  tubes  are  assumed  to  be  laid  out  in  a 60 
degree  triangular  pattern  for  aaxicua  ccapactn ess.  The 
tubes  are  located  at  the  vertices  of  an  equilateral 
triangle.  iwc  cross  sectional  areas  are  calculated.  The 
area  of  a hexagon  whose  ainor  radius  is  egual  to  cne  half 
tbe  pitch  is  deteriined  froa 

A - 6(^y  TAN  30°  (15) 


Then  the  effective  cross  sectional  area  of  the  tube  bank  is 
fcund  by  subtracting  the  clearance  area,  SA,  froa  the  cross 
sectional  area  of  the  shell  after  correcting  the  shell 
diaaeter  fer  a clearance  space  to  allow  for  free  flow 
channels  both  within  the  tube  bundle  and  arcund  the  octside 
of  the  tube  bundle.  The  effective  cross  sectional  area  of 
the  tube  bark  is  divided  by  tbe  area  of  the  hexagen,  found 
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above,  tc  determine  the  nuaber  of  tubes  in  a shell. 

lie  lieu  of  the  liquid  and  vapor  through  the  tube 
back  is  assuied  to  be  turbulent.  The  inccaing  sutcccled 
noiking  fluid  fren  the  teed  pump  lixes  witfc  the 
recirculating  fluid  returning  tren  higher  up  in  the  tube 
bank.  The  circulation  of  the  fluid  is  induced  ty  the 
boiling  process.  lhis  mixture  enters  the  tube  bank  and 
begins  the  bciling  process  Near  the  bettoa  of  the  tuba 
bark  the  beat  transfer  goes  aostly  to  bring  the  fluid  up  to 
the  boiling  teaperature.  After  the  fluid  reaches  saturated 
liguid  ccnditicns,  stable  boiling  begins.  The  vapor  bubbles 
rising  through  the  surrounding  liguid  increase  the 
turbulence  arcund  the  tubes.  As  the  liquid  continues  up 
through  the  tube  bank,  the  local  boiling  teaparatcra  is 
decreasing  due  to  the  decreasing  hydrostatic  head.  This 
decrease  in  the  boiling  teaperature  increases  the 
teiperatcre  difference  and  increases  the  heat  flew  and 
increases  tte  bubble  generaticn  rate.  Tne  vapor  generation 
■ay  becoae  sc  rapid  that  a condition  called  vapor  blanketing 
nay  cccrz  in  the  interior  of  the  tube  bank  where  the  liguid 
cannot  flew  inward  fast  enough  to  displace  the  vapor  flowing 
through  the  region.  This  is  siailar  to  burnout  in 
convective  bciling  in  that  the  tubes  dry  out  and  the  heat 
transfer  rate  decreases.  This  effect  becoaes  pronounced 
only  near  the  critical  heat  flux  according  to  Starczewski, 
in  Bef.  6. 

All  cf  the  heat  and  tas s flow  variations  in  the  tube 
bundle  are  averaged  by  performing  the  calculations  fer  a 
sirgle  tube  located  at  the  aid  height  of  the  tube  bank. 
This  tube  is  assumed  to  represent  the  average  of  all  the 
tubes  in  the  tube  bank.  The  boiling  situation  is  assuied  to 
be  saturated  pool  boiling.  Quoting  frea  Collier,  Ref  12, 
"pcol  boiling  is  defined  as  boiling  froa  a heated  surface 
subaerged  in  a large  voluae  of  stagnant  liquid".  Since  the 


fluid  vclmt  seen  ty  the  tube  is  neither  large  nor  stagnant, 
a ccrrecticn  is  reeded  for  the  effect  cf  being  in  a tube 
bundle.  Saturated  boiling  is  defined  as  boiling  at  a 
constant  temperature  and  pressure  as  the  liguid  changes  to  a 
vapor.  This  is  picbabiy  not  strictly  true  because  the 
boiling  teaperature  is  controlled  by  the  local  pressure  at 
each  tube,  and  the  fluid  is  not  reaching  a stable 
teaperature  as  it  rises  in  the  bundle.  Ihe  situation  is 
probably  closer  to  convective  boiling  since  the  fluid  is 
belling  in  a confined  channel  between  the  tubes.  In 
additional  effect  in  the  boiler  is  the  fact  that  the 
hydrostatic  bead  at  the  bottom  of  the  tute  bundle  causes  the 
boiling  tc  take  place  at  a higher  pressure  at  the  bettoa 
than  the  tciling  at  the  top  of  the  bundle  where  the  pressure 
is  equal  tc  the  boiler  exit  pressure.  This  effect  is 
accounted  fer  by  using  the  average  head  of  the  bundle  in  the 
heat  tracsfei  equations. 

Ibe  beat  transfer  analysis  of  the  toiler  begins  with 
the  basic  equation 


Q = o*  AT*,  116) 

where  U is  the  overall  heat  transfer  coefficient,  f-.  is  the 
area  tbrcugb  which  the  heat  is  transferred  and  AT*,  is  a 
properly  defined  aean  teaperature  difference.  In  this 
paper,  the  product  UA  is  treated  as  a single  quantity  called 
the  therial  conductance.  One  aay  look  at  the  equation  as  an 
analogue  cf  the  basis  electrical  equation  I * E /A  where  I is 
equivalent  tc  Q,  E is  the  saae  as  AT*,  and  UA  is  1/E.  This 
analogy  is  used  when  developing  the  expression  for  UA. 

In  ctder  to  account  for  the  effects  of  the  tube 
seeing  aany  nearby  neighbors,  the  equation  is  modified  by 
the  inclusion  of  a factor,  F,  to  account  for  the  effects  of 
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the  surrounding  tubes.  Equation  16  is  changed  to 


Q = OA  F ATm  (17) 

where  F is  defined  as  the  ratio  of  the  heat  flux  of  the  tube 
bundle  tc  the  heat  flux  for  a single  tube  in  a stagnant  fool 
boiling  situation.  According  to  Ref.  13#  a report  published 
by  workers  at  Heat  Transfer  Research,  Inc.  (HTEI)  of 
Alhambra,  California,  F is  greater  than  cne.  In  Figure  5, 
which  is  sketched  frco  a figure  in  their  report,  F ranges 
fncm  4.73  at  ATm  * 8°F  to  2.63  at  = 40°F.  Since  nc  data 

is  given  cn  the  conditions  of  the  test.,  F is  assumed  to  be 
fixed  at  three. 

Ir  crder  to  find  AT„,,  the  heat  exchanger 
effectiveness,  6 , is  determined.  The  effectiveness  is 
defined  as  tte  ratio  of  the  actual  heat  transfer  rate  to  the 
uaximum  possible  neat  transfer  rate  as  shown  by 

s C6  (Th,in  ~ T\out  ) 
c*v»iw  (rh,in  " Tc,in  ) 

where 

C = («  (19) 

ctwih  = °f  Ch  or  Cc  (20) 

In  the  fccilei,  Th  is  the  temperature  of  the  entering 
seawater,  ard  TC|in  is  the  saturation  temperature  cf  the 

working  fluid  at  the  pressure  at  the  mid  height  of  the 
boiler.  Reference  14,  by  Kays  and  London,  contains 

equations  for  finding  the  ef f ecti veness  for  various  flow 
configurations.  In  addition  tc  finding  Cmin , cmax  must  be 


FIGURE  5.  COMPARISON  OF  A SINGLE  TUBE 


kncwn  in  crder  to  decide  which  equation  applies.  For  the 
case  of  tciling,  Cmax  is  infinite;  therefore,  regardless  of 
the  flew  pattern 

i = 1 - BXP(-NIU)  (21) 

where  NIC,  the  number  of  heat  transfer  units  cf  a heat 
exchanger,  is  defined  as 

NTO  = UA/Cmin  (22) 

and 


Cmin  = 


(22) 


AT*,is  found  fren 


£ 


UA  AT, 


•mm 


AT, 


(24) 


and  the  temperature  of  the  sea  water  at  the  outlet  cf  the 
tube  is  deteraiued  frea 

£ * ~ ^tyouf  (2C) 

where 

AT0  = lh<in-  Tc,in  <26) 

lhe  quantity  UA  is  deterained  next.  Figtre  6 
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FIGURE  6.  THERMAL  NETWORK  FOR  A 
SINGLE  TUBE 


39 


illustrates  the  method  used.  If  UA  is  pictured  as  the 
reciprocal  cl  resistance  then  the  flow  cf  energy  from  one 
fluid  tc  tie  ether  can  be  pictured  as  a series  electrical 
circuit  and  U A is  determined  from 

OA  = 1/  (R 1 ♦ R2  83  ♦ R4)  \21) 

Reference  15,  by  Holman,  shows  how  each  resistance  can  be 
feund.  El  is  the  thermal  resistance  due  tc  convecticn  on 
the  inside  cf  the  tubes,  the  sea  water  side. 


r" 


i 


R 1 


1 


7T  h 


sw 


L 


R2  is  the  thermal  resistance  of  the  tube  wall. 


(28) 


ln(d0/d,  ) 
2 7T  L 


(29) 


R3  is  tie  thermal  resistance  attributed  tc  boiling  cn  the 
outside  cf  tie  tubes. 


R3 


1 


(30) 


R4  is  the  thermal  resistance  due  to  any  other  resistances, 
such  as  feuling  and  corrosion  deposits. 


Tie  beat  transfer  coefficient  on  the  seawater  side, 
h5w , is  determined  from  the  dimensionless  group  called  the 
Nusselt  number  defined  as  follows  for  flow  inside  of  tubes. 
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mmw  r,; — 


Nu 


hsw  d; 


(31) 


In  cider  tc  find  the  Nusselt  number,  three  distinct  regions 
cf  flow  are  allowed  for  in  the  problem:  laiinar, 

transiticr,  and  turbulent  flow.  The  division  points  are  in 
terms  of  valres  of  the  Reynolds  number.  Re  = ud p/p  , 

laminar:  0 < Re  < 2,(JC0 

transition:  2,000  < Re  < 10,000 

turbulent:  10,000  < Re  < infinity 

In  the  laiinar  region  the  Siedec-Tate  correlation  is  used. 

1/3 

NUL  = 1.86  (Be  Pr  dj /L  ) (32) 

In  the  turbulent  region  the  Littus-Boelter  correlation  is 
used. 


C.8  0.3 

KuTu  = 0.023  Re  Pr  (33) 

In  the  transition  region  the  flew  is  unstable,  but  froi  data 
in  Ref.  14  it  appears  that  a fairly  smooth  curve  drawn  from 
the  laminar  curve  to  the  turbulent  curve  adaguately 
represents  tbe  value  cf  the  Nusselt  number  in  the  transition 
regicn.  Ibis  is  done  by  assuming  the  form 

■«Tr  = y mul  ♦ (1  - y)  nuTm  (34) 


where 


3 2 

y » a He  ♦ t He  ♦ c Se  ♦ d 


(35) 
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The  coefficients  in  the  polynomial  are  found  ty  reguiring 


fv 

\ 


NuTr  = NuL  at  Re  = 2,000 
KuTr  = NuTu  at  Re  = 10,000 

ana  that  the  first  derivatives  with  respect  to  the  Reynolds 
muter  natch  at  Re  egual  to  2,000  and  10,000. 

The  I. I.  Mostinski  correlation,  from  Ref.  16,  for 
sirgle  hcrizcrtal  tube  pool  boiling  is  used  to  find  the  heat 
transfer  coefficient  on  the  outside  of  the  tube.  HTEI's 
extensive  experimentation  led  Palen,  Yarden,  and  Tatcrek,  in 
Ref.  13,  tc  conclude  chat  Mostinski’s  eguaticn  produced  the 
"»cst  consistent  results".  His  correlation  is 


C . 6 9 0.7  0.17  1.2  1C 

hwj  = C.CC658  pc  g (1.8  r ♦ 4 r ♦ 10  r ) (36) 


where 


r = P/Pc 

Pc  * critical  pressure  of  the  working  fluid  (psia) 
q = heat  flux  (BTO/hr-f t*) 


A final  heat  flux  relation  needed  is  the  relation  fer  the 
maximum  heat  flux  allowed.  Mostinski, s correlation  for  a 
sirgle  tube  is  used. 


0.35  0.9 

<3m<nf  * 303  pc  r (1  - r)  (37) 

where  the  urits  are  the  sane  as  in  eguaticn  (36).  Froa 
HTEI’s  data,  under  unknown  conditions,  gwaJf  for  the  tundle 
appears  tc  he  one  third  of  the  g^a)(  for  a single  tube. 
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lhe  fewer  required  tc  pump  the  salt  water  tbreugh 
the  tubes  depends  cn  the  pressure  drop  across  the  tube  bank. 
There  are  three  causes  of  the  pressure  drop.  The  flow 
experiences  a loss  as  the  water  enters  into  and  exits  from 
the  tubes  in  addition  to  the  frictional  loss  along  the  tube. 
Free  Ref.  14,  the  pressure  drep  fer  a liquid  is  determined 
frci 


2 4 f L 

Ap  * 0.5  pgw  Uj*,  ( — ♦ Kc  - Ke)  (38) 

di 

In  erdez  tc  determine  f,  Ke  and  Kc,  the  flow  regime  must  be 
specified.  The  regimes  are  a little  different  than  those  for 
the  Nusselt  cumber,. 

laminar:  0 < Re  S 2,0C0 

transition:  2,000  < Re  £ 5,000 

turbulent:  5,000  < Re  < infinity 

The  Fanning  friction  factor,  f,  in  the  laminar  regime  is 

f « 16/Re  (39) 

In  the  turbulent  regime,  the  Fanning  fricticn  factcr  is 
found  fret  the  Bias  's  eguation  in  the  range  5,00C  < Re  < 
ICC, 000 

0.25 

f - 0.079/le  (40) 

which  is  a close  approximation  to  the  Karman  - Hikcradse 
eqcaticn , 
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-0.5  0.5 

<4  f)  = -0.6  ♦ 0.87  In  (R € (4  f)  ) 


<4  1 ) 


which  is  valid  throughout  the  turbulent  region  (Re  > 5,000), 
assuming  snccth  walled  tubes.  Either  eguaticn  (40)  cr  (41) 
■ay  be  used  with  siiilar  results  but  the  Elasius  is  tetter 
suited  tc  this  analysis  because  the  Karman-Nikuradse 
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The  entrance  effect,  Kc,  and  the  exit  effect,  Ke, 
depend  on  the  Beynclds  number,  the  free  flow  to  frontal  area 
ratio,  4>  , and,  in  the  laminar  region,  on  the  tube  length. 
Fcr  the  toiler,  $ is  determined  by 


4> 


d;  n 


Ds  - 4 SA 


(42) 


Figure  S,  page  7-13,  in  Ref.  17  is  used  tc  find  Ke  and  Kc. 
The  effect  cf  Reynolds  number  and  L/d  ratio  was  eliminated 
from  tbe  relationships  in  the  laminar  and  the  turtulent 
regimes  cf  the  flow,  tut  the  Reynolds  number  was  included  in 
the  trarsiticn  region  because  of  the  large  variation  of  Kc 
with  Reynolds  numter.  The  following  eguaticns  are  used  to 
find  the  entrance  effect  in  the  various  flow  regimes. 

laminar:  (using  the  curve  for  4 (L/d) Re  - 0.1) 


Kcl  ■ 1.  - 0.4  6 


(43) 
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Turbulent:  (using  the  curve  for  he  = 5,0GC) 


KcTu  = 0.52  - 0.4  <} 


(44) 


Iransiticr : 


KCjy  = A KcL  ♦ (1  - A ) KCju 


(45) 


where 


A * e Be  ♦ f (46) 

The  ccelficients  e and  f are  determined  by  requirirg  the 
transition  Kc  to  egnal  the  laminar  Kc  at  Be  equal  to  2,000 
and  the  turbulent  Kc  at  Be  equal  to  5(000.  The  exit  effect 
is  the  saae  fer  all  flow  regimes,  and  is  found  frci  the 
fcllcwing  equation  that  approximates  the  4 (L/d) Be  =0.05 
curve  in  Eef.  17. 

2 3 

Ke  = 1.CC93  - 2.5178  ^ ♦ 1.1613$  - 0.  17677  $ (47) 


2.  SSBdfifliSI 


Its  design  of  the  condenser  is  almost  the  saie  as 
the  design  cl  the  boiler.  It  is  a fixed  tube  sheet, 
tute-in-stell  beat  exchanger.  The  condenser  does  net  reed  a 
separation  velum*  mo  the  entire  shell  is  filled  with  tubes. 
In  erder  to  allow  the  entering  vapor  to  reach  all  cf  the 
tubes,  channels  through  the  tube  bank  are  needed.  The 
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allowance  fo£  this  is  accomplished  by  decreasing  the  shell 
diaieter  by  a fixed  factor  when  calculating  the  number  of 
tubes.  The  calculation  of  the  working  fluid  side  heat 
transfer  coefficient,  h , is  different . The  equation 
acccunts  for  the  condensate  from  one  tube  falling  cn  those 
below  it  ard  for  the  subcooling  of  the  condensate.  M.M. 
Chen's  equation,  as  developed  in  Ref.  18  and  shown  in  a wore 
siiple  ferr  in  Ref.  17,  is  used. 


h*f  = C .7  26  fl+0. \ (n-  1)1  ll-fi  . } * Hffl 

f L \«f0  / 'JL  n de/4  AI  J 


0.25 


(48) 


where 


k = BlU/f t2-hr-°F 


- 1Sot  - T*w  <“*> 

3 

H!  * H,  ♦ c„AT  (50) 

*9  8 p 

n = number  of  tubes  in  a column 

is  the  enttalpy  corrected  for  subccoling  of  the 
ccrdensate.  Equation  (48)  replaces  equation  (36)  ir.  the 
analysis  of  the  condenser.  In  the  circular  shaped  tube 
bark,  n is  found  by  determining  the  number  of  tubes  in  a 
column  of  average  fceight. 


Ds 

n = 9.425  1 (51) 

Pi 


The  definition  of  ATC  becomes 
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* 1*ot  “ 


in  contrast  tc  equation  (26)  above, 
teaperatuie  is  found  froa 


Ihe  seawater  cutlet 


£ = Tsw/out  ~ Tsw,»n 
AT„ 


which  replaces  equation  (25) . 

In  order  tc  find  the  heat  transfer  ccefficient,  the 
teaperatore  cf  the  outside  tube  surface  is  needed  in 
equation  (49);  therefore,  the  theraal  network  is  solved  for 
tfce  outside  tall  teaperature  and  for  the  heat  flow  rate. 


04*  (T„  - T^) 


where 


UA'  * 1/ (F 1 ♦ E2  ♦ R4)  (55) 

^bulk  * aTera9e  °f  the  inlet  and 
the  outlet  seawater  teap. 

There  is  nc  factor,  F,  in  the  heat  transfer  eqnaticn  for 
ccrdensation  since  the  effect  of  neighboring  tubes  is 
accounted  fcr  in  equation  (47) . 

3* 


Eoth 


Here  are  three  puaps  in  the  aodel.  Eoth 
circulation  puaps  for  the  boiler  and  the  condenser  are 
assuaed  to  be  propeller  puaps  (axial  flow  puaps).  The  feed 


^7 


pumps  are  assumed  tc  be  centrifugal  pumps.  The  differences 
in  the  characteristics  of  the  tec  types  of  pump  show  up  cnly 
in  the  calculation  of  capital  costs  of  the  pumps  ir  the 
model . 


lie  choice  of  the  type  of  pump  to  apply  is  cade  on 
the  basis  cf  the  flow  rate,  GEM,  and  on  the  head  the  pump  is 
working  against.  Cualitatively,  the  propeller  pump  is  used 
where  the  flew  rates  are  high  and  the  head  is  lew.  The 
centrifugal  pumps  apply  where  the  flow  rate  is  comparatively 
lew  and  the  head  pressure  is  high.  The  qualitative  picture 
is  placed  in  tetter  perspective  through  the  calculation  of 
the  specific  speed  from 


Nsg  = 


0.5 

fi  PH  GEM 
0.75 

Head 


(56) 


The  centrifugal  pump  is  used  when  500  < Nsg  < 7,500  and  the 
propeller  puip  applies  when  7,500  < Nsg  < 15,000,  acccrding 
to  Eef.  IS.  The  choice  of  the  type  of  pump  is  checked  after 
the  optimization  process  has  been  completed  by  sclving 
eguation  (55) . If  the  value  of  the  specific  speed  is 
outside  cf  the  range,  500  < Nsq  < 15,000,  the  maximum 
efficiency  cf  the  pump  drops  off  rapidly;  therefore,  few 
pumps  are  designed  tc  operate  outside  cf  this  range. 

Since  the  pumps  are  "black  boxes",  only  the  power 
consumed  by  the  pumps  is  calculated.  The  power  required  to 
drive  the  pump  is  found  from 


M Ap 


(57) 


c. 


COST  FFAREhOEK 


As  discussed  at  the  start  of  this  report,  the  capital 
costs  of  the  major  components  are  the  only  costs  to  be 
considered  ir  this  mcdel.  The  boiler  and  the  condenser  are 
both  tute-in-snell  heat  exchangers  assuied  to  be  cf  the 
fixed  tube  sheet  type.  The  only  difference  being  that  the 
shell  of  the  boiler  is  partly  filled  Kith  tubes.  Eoth 
seawater  circulation  pumps  are  assumed  to  be  propeller  type 
puips  and  made  cf  the  same  material.  The  size  and  number 
per  heat  exchanger  shell  may  be  different  fcr  the  toiler  and 
the  ccnaensei.  The  number  of  working  fluid  feed  pumps, 
which  are  the  centrifugal  type  pumps,  are  assumed  tc  egual 
the  number  cf  toiler  shells. 

Ihe  ccst  estimating  relationships  fcr  the  heat 
exchangers  and  the  centrifugal  pumps  are  taken  from  Pef.  20, 
by  K.H.  Guthrie.  The  eguation  fcr  the  cost  of  the  heat 
exchangers  is 


*h<  * 116.3  (Fd*  F„)  F^  1^  (At)°*63  150 

where  the  constants  are  calculated  from  a leg-log  plot,  of 
cost  versts  total  tube  surface  area  and 

Fd  s 0.6  for  a fixed  tube  sheet 

Fp  = C.C  for  pressure  correction  (p  < 150  psi) 

F^  2 material  factor  from  Table  1 
Ihe  2 ccst  index 
A<  * de  L N 

fcr  a relevant  range  of  100  < A < 10,000  ft*. 

The  ccst  eguation  for  the  centrifugal  pumps  and  the 
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associated  meters  is 


0.721 

*fp  = ^C13  Fo  CH  (59) 

where  the  constants  are  found  from  a straight  line 

approximation  cf  the  log-log  plot  cf  ccst  versus  C for  the 

section  cf  the  curve  from  C = 30,000  to  C = 300,000  fcr  an 
electric  motor  driven  pump  and 

F0  = 1.0  fcr  suction  temperatures  < 250°F  (121°C) 

Ipp  = ccst  index 

Fm  = material  factor  from  Table  2 
C j_|  = predict  of  the  flow  rate  in  GFM 

and  the  pressure  differential  in  psi 

fcr  a relevant  range  of  30,000  < CH  < 300,000.  The  cost 

data  in  Guthrie's  article  is  given  for  a time  fcase  of 

■id-1968. 

Feference  21,  also  ny  Gutnrie,  gives  cost  data  for 
propeller  pumps  and  centrifugal  pumps.  The  eguation  fcr  the 
propeller  pumps  with  the  motors  is 


0.783 

Scp  * 5.767  Fp  lcp  Flow  (60) 

where  the  constants  are  found  from  the  log-leg  plot  cf  cost 
versus  flew  rate  for  an  electric  motor  driven  pump  and 

Fp  = 1.0  fcr  a suction  pressure  < 150  pei 

F^  = material  factor  from  Table  3 

Itp  * cost  index 

Flow  = flew  rate  in  GPn 

for  a relevant  range  of  1,000  < Flow  < 100,000  GEK.  The 
time  base  fcr  the  propeller  pump  data  is  fcr  the  erd  of 


Table  1 


flJTERJAI  FACTORS  FCB  TH £ HEAT  EXCHANGERS 


Material 


Shell 

Tube 

F f 

Carter  Steel 

Carbon  Steel 

1.00 

Carton  Steel 

Brass 

1.52 

Carton  Steel 

Stainless  Steel 

3.  52 

Stairless  steel 

Stainless  Steel 

4.50 

Carton  Steel 

Monel 

3.  75 

bene  1 

Monel 

4.95 

Carbon  Steel 

Titaniua 

11.10 

Titariun 

Titanium 

16. 6C 

Table  2 

FATEEIAL  FACTORS  FOR  THE  CENTBIFUGA1  FEED  E UFF 


Batetial 

fJS 

Cast  Iron 

1.00 

Bronze 

1.28 

Cast  Steel 

1.32 

Stainless  Steel 

1.93 

Monel 

3.23 

Hastelloy  C 

2.89 

Titaniua 

8.  98 
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Table  3 

MATERIAL  FACTORS  FOR  PROPELLER  PUMPS 
pumps 


Material 

E n 

Cast 

Iron 

1 .00 

Cast 

Steel 

1.28 

Stai 

nless  Steel 

1 .64 

52 


' - - 1 — 


1570. 

lie 

costs 

arri 

ved  at 

froa  egu 

atior 

£ (58) 

9 

(59)  , 

and 

<6C) 

are 

in 

dollar 

s per 

coaponent,  for 

exam 

pie,  dol 

lars 

per 

heat 

ex 

chai 

nger  s 

bell. 

The 

purpose 

of 

the  cost 

ince: 

x is 

allow 

al 

1 the  egui 

paent 

to  be 

costed  o 

n a c 

c it  oc  n 

ba 

sis . 

Tne 

ties 

di 

f f€] 

ter.ee 

is  a 

ssuaed 

to  be  sa 

all  e 

ncugh 

so 

that 

the 

cost 

ind 

ex  : 

is  set 

to  1 . 

00  for 

all  the 

compo 

nents. 

E.  ASSEPEI1  CP  THE  MATHEMATICAL  MODEL 

The  engineering  and  the  ccst  relationships  are  linked 
together  in  the  acdel.  The  cost  equation  is  the  function  tc 
be  ainitized  and  the  implicit  and  explicit  constraints  are 
the  engireerirg  eguations  that  restrict  the  feasible  region 
of  the  (ictlea.  The  explicit  constraints  restrict  the 
variables,  the  vector  X,  to  positive  values.  The  rest  of 
the  BOdel  is  fcriulated  in  the  following  fcri 

■ iriaize  f (X) 

subject  to:  g^(X)  > 0,  i=  1,2 , . . . . ,a 

h (X)  = 0,  i*a*1 , . . . . ,a+n 
i 

where  tie  definitions  of  the  elenents  of  X are  shown  in 
Table  4.  Several  of  the  variables  Bust  have  only  integer 
values  if  the  systea  were  actually  built,  but  for  the 
objectives  cf  this  acalysis  they  are  left  as  continuous  real 
variables.  Per  exanple,  the  nuaber  of  shells  aay  have  an 
integer  value  in  the  actual  systea  but  this  reguireaert  can 
be  accoaccated  by  solving  the  systea  for  a real  value  then 
integerizing  the  nuaber  of  shells  and  solve  the  prcblei  with 
the  next  higter  anc  the  next  lewer  integer  nuaber  of  shells. 
The  lower  cf  the  two  solutions  would  be  chosen.  lhe  final 
resolution  cr  the  correct  nuaber  of  shells  would 
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Table  4 

DEFINITION  CF  THE  ELEMENTS  OF  THE  VECTOR  X 

x (1 ) = velocity  of  seawater  in  the  fcciler 

x(2)  = " " " " " condenser 

x (3)  = cutlet  pressure  in  the  boiler 

x(4)  = inlet  " " " condenser 

x(5)  = tube  length  in  the  boiler 

x (6)  = " " " " condenser 

x (7)  = lass  flow  rate  of  the  wcrking  fluid 

x(6)  = diaaeter  of  the  bciler  shell 

x(S)  = " " " condenser  shell 

x(1C)  = outside  dia.  of  a tube  in  the  boiler 

x(11)  * " " ii  ii  n ii  ii  condenser 

x(12)  = nunber  of  bciler  shells 

x(13)  = *'  •'  condenser  shells 
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prcbatly  depend  on  other  than  thermody namic  and  cost 
consider aticcs , such  as  a requirement  that  cne  extia  shell 
be  included  to  allow  for  one  shell  being  dcwr  for 
maintenance.  Obviously,  the  number  of  tubes  in  a shell  must 
be  an  irtecer  value  but  cne  extra  tube  cut  of  several 
thousand  is  a minor  matter.  In  addition,  there  may  be  a 
requirement  fcr  extra  tubes  to  allow  for  a percentage  tc  be 
plugged  wten  they  develop  leaks  so  that  the  entire  shell  is 
net  put  cut  cf  ccmmission  by  a minor  leak  in  one  tube. 

1.  Objective  Function  and  the  Constraints 

The  objective  function  is  the  sum  of  the  costs  of 
the  boiler  arc  the  condenser,  and  the  pumps.  The  model 
assumes  that  there  is  one  circulating  pump  and  one  feed  pump 
for  each  shell  in  the  bcilez  and  that  there  is  one 
circulating  pump  for  each  shell  in  the  condenser.  Since 
there  is  nc  cost  analysis  for  the  turbine,  nc  assumption  is 
needed  about  the  number  of  turbines  in  the  system. 

There  are  nineteen  constraints  in  the  present 
pregram,  cne  cf  which  is  optional.  The  first  three  concern 
the  saturation  pressures  in  the  boiler  and  the  condenser. 
The  next  twelve  place  restrictions  on  various  dimensions  of 
the  beiler  anc  the  condenser.  The  last  three  reguire  that 
seme  intermediate  calculations  stay  positive. 

Csing  the  objective  function  as  a constraint  has 
proved  advantageous  in  cases  where  the  problem  was 
ncn-ccnvex  or  where  the  program  had  difficulty  in  finding  a 
solution.  It  is  used  by  setting  a boundary  cn  the  objective 
function  sc  that  the  constraint  is  cither  infeasible  or 
slightly  feasible  at  the  starting  point.  The  proper  choice 
must  be  feund  by  trial  and  error,  one  caution  to  note  is 
when  equality  constraints  are  involved.  If  the  starting 
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pcint  1 s rtai  the  optimum,  the  program  tries  to  move  away 
from  the  boundary  and  the  equality  constraints  may  not  be 
satisfied  when  the  program  finishes. 

lie  need  for  the  ether  eighteen  inequality 
constraints  is  net  obvious  when  the  problem  is  first 
formulated.  The  thermodynamics  of  the  problem  require,  for 
exaiple,  that  the  propane  saturation  temperature  be  less 
than  seawater  temperature  in  order  for  the  bciler  tc  absorb 
energy,  bet  the  program  has  no  way  of  knowing  that,  unless 
it  it  specified.  If  the  program  can  impreve  the  cost  and 
satisfy  the  ether  constraints  by  reversing  the  relationship 
it  will  do  so.  This  inversion  may  take  place  only  for  a 
shert  time  doting  the  solution  of  the  problem,  but,  if  the 
equations  fail  under  this  condition,  the  constraint  on  the 
temperatures  in  the  bciler  is  necessary.  best  of  the  other 
constraints  ace  needed  for  the  same  reason. 

Seme  constraints  are  net  added  to  the  problem  since 
they  are  net  bindinq  in  the  solution.  There  are  twe  such 
ccnstraints.  The  first  is  the  segmental  area  reguired  to 
allow  separation  of  the  liguid  droplets  from  the  vaper.  The 
segmental  area  set  aside  by  the  4055  rule-of-thumb  is  more 
thac  that  required  by  the  consideration  of  the  physical 
preperties  cf  prepane.  The  ether  constraint  is  the 
restriction  cn  the  maximum  heat  flux  by  Nostinski's 
correlation,  equation  (37).  The  constraints  dc  not  become 
binding  in  the  solution. 

The  three  equality  ccnstraints  contain  the 
calculations  cf  the  properties  of  the  various  components 
modeled  in  this  analysis.  The  first  one  concerns  the 
bciler.  It  states  that  the  actual  energy  absorption  rate  of 
the  bciler  shells  must  be  equal  tc  the  energy  absorption 
rate  reguired  by  the  thermodynamic  cycle.  Similarly,  the 
second  eguality  constraint  states  that  the  beat  rejected  by 
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the  condenser  must  equal  that  reguired  by  the  thermod ynamic 
cycle  for  the  condensation  of  the  vapor  from  state  point  5 
back  tc  state  pcint  1 as  shown  in  Fig.  3.  The  third 
constraint  calculates  the  pcwer  required  by  the  various 
pumps  and  the  pcwer  extracted  from  the  turbine  to  run  the 
system's  pumps  and  to  produce  electrical  power  for 
consumption.  Ihis  ccnstraint  requires  that  the  power  cut  cf 
the  turbine  equal  the  power  used  by  the  pumps  anc  the 
generator.  The  net  useful  power  cut  of  the  system  is  a 
design  objective  and  is  specified  as  25MH  herein. 

2.  2 be  Optimization  Method 

Sirce  the  arrival  of  the  electronic  computer,  a 
great  many  nonlinear  programming  techniques  have  teen 
developed  from  early  theories  and  more  current  research. 
The  computer  program,  called  the  Sequential  Unconstrained 
Minimization  Technique  (SUHT),  developed  by  Hylander, 
Helens,  and  McCormick  in  Be f.  22  is  used  to  optimize  the 
mcdel.  Their  pregram  implemented  much  cf  the  theory 
contained  in  fief.  23,  by  Fiacco  and  McCormick,  about 
ncnlinear  programming  using  unconstrained  minimization 
techniques.  SUMT-fersion  4 is  chcsen  because  of  its  ability 
to  solve  a wide  variety  of  problems.  The  theoretical 
requirements  for  SUHT  to  be  guaranteed  to  find  a local 
minimum  is  that  there  exists  seme  point  which  satisfies  the 
inequality  ccnstraints  and  that  there  not  be  a local  minimum 
at  pcints  where  X goes  to  infinity.  If  certain  convexity 
ccrditicrs  are  satisfied,  then  the  local  minimum  found  is  a 
glcbal  minimum. 

£ GMT  uses  a penalty  functicn  tc  transform  the 
ccnstraircd  optimization  problem 
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f (X) 


subject  tc:  g (X)  > 0,  i=1,2,....,m 

i 


h (X)  = 0,  i=  *♦  1 , . . . . , men 

i 


intc  a sequence  of  unconstrained  problems  cf  the  form 


F(X,r)  = f (X)  - r I In  g (X)  ♦ 1/r."in  (h  (X) ) 2 (61) 

1=1  l 1=  me  1 l 

where  E|X,r)  becomes  the  function  to  be  minimized.  The 
inequality  constraints,  g (X)  , cauje  a large  penalty  to  occur 
as  a boundary  is  approached.  The  equality  constraints  add  a 
large  penalty  whenever  h(X)  departs  from  zero  in  either 
direction.  The  size  of  the  penalties  due  tc  g (X)  and  h ( X) 
is  controlled  by  the  arbitrary  parameter  *'r".  The  algorithm 
solves  the  problem  by  solving  a series  of  sukproblems.  For 
r1  >r^>r3  . . .>0,  SOHT  minimizes  E(X(r„),rn)  in  each  sutprotlem 
in  order  to  find  the  minimum  of  f (X)  as  "r"  tends  to  zero. 
When  a minimum  of  a subproblem  is  found,  "r"  is  reduced  by  a 
factor  assigned  by  the  user  and  the  new  subprcblem  is 
solved.  The  process  continues  until  the  convergence 
criteria  chcsen  by  the  user  is  satisfied.  Figure  7 
illustrates  schematically  the  results  cf  the  sequerce  of 
sutprcblems  and  the  qualitative  shape  cf  the  penalty 
furction.  Three  methods  of  minimization  are  included  in  the 
program.  The  user  may  choose  either  the  generalized 
Newton-Faphscn,  steepest  descent,  or  McCormick  modification 
of  the  Fletcber-Po well  method. 

There  are  several  other  options  that  contrcl  the 
operation  of  SHUT.  The  user  may  cheese  the  initial 
penalty,  r,  cr,  if  the  problem  satisfies  certain  conditions, 
the  user  may  elect  to  have  SUMT  compute  the  initial  "r". 
There  are  three  choices  for  the  convergence  criterion  cf  the 


58 


X 

FIGURE  7.  PROBLEM  SOLUTION  SEQUENCE 
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problem.  The  criterion  used  depends  on  the  inequality 
penalty  ter*  being  close  enough  to  zero.  Another  option  is 
the  choice  of  the  calculation  of  the  subprcblem  stepping 
criterion  . Two  of  the  choices  compute  the  inverse  cf  the 
seccnd  partial  matrix  and  the  third  choice  computes  only  the 
gradient  cf  the  penalty  function  ( VxP(X,r)).  One  cf  the 
acst  useful  cptions  is  extrapolation.  After  two  or  three 
subprofcleis  have  been  solved  SUHT  is  able  tc  make  a first  or 
seccnd  crcer  extrapolation  to  the  starting  pcint  for  the 
next  suhprcbleo.  This  can  speed  up  the  solution  cf  the 
problem  meascreably. 

Ike  characteristics  of  the  problem  limit  the  choices 
for  each  cpticn.  Ihe  penalty  factor  couldn't  be  set  by  SUflT 
due  to  the  presence  of  equality  constraints  and  because  the 
iritial  1 is  usually  not  close  to  any  boundary.  The  stepping 
criterion  fer  the  problem  is  that  the  inequality  penalty 
ten  beccme  less  than  some  value,  usually  10~7. 
Extrapolation  is  used,  but  it  makes  little  difference  in  the 
answer  or  tc  the  solution  time  whether  a twe  pcint  cr  three 
point  estiiate  is  made.  The  subprctlem  convergence 
criterion  is  the  magnitude  of  the  gradient  of  P(X,r) 
becoming  less  than  some  small  number. 

The  choice  cf  starting  penalty,  r,  affects 
significantly  the  solution  time  and  even  whether  a solution 
is  found  it  seme  problems  according  to  Pef.  7.  If  the 
initial  "r"  is  toe  small,  the  subproblem  optimum  may  be  the 
problem's  optimum  which  may  be  hard  to  find  cr  require  more 
calculations  than  if  a much  larger  "r"  is  chcsen.  If  "r"  is 
toe  large,  tbe  first  few  optima  lie  near  the  center  of  the 
feasible  region  and  are  unrelated  to  the  optimum  cf  the 
objective  function. 

Scaling  of  the  variables,  constraints,  and  the 
objective  function  is  necessary.  The  variables  must  be 
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within  twc  cr  three  orders  of  magnitude  of  each  other,  for 
eiaaple,  the  working  fluid  flow  rate  is  of  the  order  of  107 
whereas  the  tube  diaaeters  are  of  the  order  10®.  If  the 
finite  difference  interval  is  10~4,  the  partial  of  the 
eguaticns  with  respect  to  the  aass  flow  rate  nay  be  lost  due 
tc  the  roundoff  error  providing  the  only  significant 
figures.  hith  the  constraints  and  the  objective  function, 
unscalinc  may  result  in  one  tera  in  the  penalty  function 
being  dcainant,  and  the  algoritha  finds  that  nost  cf  the 
decrease  in  the  penalty  function  takes  place  in  the  dcainate 
ter*.  The  result  is  that  SUHT  works  hard  tc  decrease  that 
tern  and  does  not  see  the  other  teras  resulting  in  a 
meaningless  solution. 

In  Baking  first  and  second  order  estimates  cf  the 
solution,  i: C H T aay  generate  negative  nutters  in  the  X 
vector.  Wfcen  it  evaluates  the  constraints  and  objective 
function  for  this  point,  the  eguations  fail  and  aay  cause 
the  program  tc  fail  depending  on  the  severity  of  the  failure 
and  how  many  tiaes  the  error  has  cccuied.  This  problem  is 
ciicuavented  by  testing  the  vector  upon  entry  in  tc  the 
subroutine  where  the  objective  function  and  the  constraints 
are  evaluated  and  setting  the  constraint  being  evaluated 
equal  to  scae  negative  constant. 

The  vain  prograa  of  SUHT  is  altered  so  that  the 
prcblea  can  be  restarted  at  the  current  scluticn  without 
reading  in  a new  problem.  This  is  necessary  because  the 
pregraa  would  give  all  the  indications  of  having  feund  a 
solution,  but  if  the  problem  is  restared  it  could  oftec  find 
a lew  value  for  the  objective  function.  The  new  solution  is 
compared  to  the  previous  solution  and  if  the  two  solutions 
are  not  close  enough  the  problem  is  restarted  after  aaking 
scae  changes  to  particular  parameters,  opticus,  or 
ccnstrairts  depending  on  the  objective  the  user  has  in  aind. 


HI.  RESULTS  AjjD  CONCLUSIONS 


One  cf  the  major  tasks  in  the  development  of  a nonlinear 
model  is  the  linking  of  the  model  with  the  optimization 
program  in  order  tc  analyze  the  model.  The  selection  cf  the 
proper  nonlinear  programming  technique  is  essential  tc  the 
solution.  Several  techniques  require  that  the  objective 
furction  and  the  constraints  be  in  a particular  form.  For 
example.  Geometric  Programming  developed  by  Euffin,  Peterson 
and  Zener,  in  Ref.  24,  requires  that  the  equations  be 
pcsyncmials  which  is  not  true  of  the  model  developed  it  this 
paper.  SUNT  was  chosen  because  cf  its  wide  applications  to 
prctlems  cf  any  type.  It  only  requires  that  there  exist 
some  feasible  point  and  that  local  minima  not  exist  at 
infinity.  It  has  successfully  solved  many  ncnconvex 
problems . 

Two  versions  of  the  analytical  model  are  presented.  The 
first  is  called  the  iterative  version  because  the 
calculation  cf  the  heat  absorption  rate  and  heat  rejection 
rate  require  iterating  on  some  quantity  to  find  them. 
Mostinski's  correlation,  equation  (36),  contains  the  heat 
flux  which  is  an  unknown.  It  is  not  possible  to  solve  the 
system  cf  equations  explicitly  for  the  heat  flux;  therefore, 
a successive  approximation  technique  is  used.  The 
calculation  cf  the  heat  rejected  in  the  condenser  involves 
the  same  process  with  the  temperature  of  the  outside  tube 
wall  as  the  iterating  variable.  SUNT  is  urable  to  arrive  at 
a consistent  solution  in  the  time  allowed,  up  to  90  minutes 
of  computer  time.  The  solutions  are  guite  different  for 
each  try  with  nc  pattern  evident  about  the  direction  the 
scluticn  might  lie. 
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There  are  several  possible  explanations  for  this 
behavior.  Perhaps  the  behavior  of  this  version  cf  the 
problem  is  such  that  SUNT  is  not  the  best  choice  fcr  the 
optimizing  algcritha  or  possibly  the  problem  cannct  be 
solved  at  all  in  this  fora.  Another  source  of  problems 
could  be  the  nuaerical  differentiation  procedure.  Nuaerical 
di  f ferentiaticn  is  subject  to  severe  truncation  and 
rcund-cff  error.  The  subproblem  always  terminates  when  SUNT 
is  unable  to  reduce  the  penalty  function  during  the  next 
iteration,  frca  one  estimate  of  X to  the  next  estimate. 
Beieaber  the  ccnvergence  criterion,  never  satisfied  for  this 
version  cf  the  problem,  is  the  magnitude  of  VXE  (X,r) 
becoming  less  than  soae  small  number.  No  conclusive  results 
could  be  obtained  with  this  version  so  the  second  version  is 
csed  fcr  the  present  analysis. 

In  this  second  version  the  heat  transfer  coefficients  on 
the  Propane  side  of  the  tube  are  fixed.  This  enables  the 
eguality  constraints  for  the  boiler  and  the  condenser  to  be 
calculated  to  the  full  accuracy  of  the  computer  without 
increasing  the  time  required.  The  correlations,  eguation 
(36)  and  (48),  are  used  for  comparison  with  the  fixed  value. 
If  the  calculated  value  and  the  fixed  value  are  within  20  tc 
2SJ  of  cne  another,  the  solution  would  be  considered  to  be 
gccd  since  this  is  the  tolerance  claiaed  for  the 
ccrrelaticns . 

SLUT  seeied  to  be  able  to  solve  the  new  prcblea  because 
a series  cf  seven  solutions  have  been  produced  where  the 
values  of  the  elements  of  X and  cf  the  objective  function 
vary  nc  aore  than  3.7%  and  1.2%,  respectively,  about  the 
acan  when  different  starting  points  and  different  SUNT 
ccntrcl  paraaeters  are  selected.  See  Table  6 for  a suaaary 
of  the  solutions  obtained  froa  the  assumed  data,  as  shewn  in 
Table  5.  The  design  implications  of  these  results  ate  guite 
reasonable.  The  range  shows  the  spread  in  the  scluticns. 
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Table  5 

PROBLEM  INPUT  DATA 


Fouling  Thermal  Resistance 

0.0  hr-°5/BTU 

Eciler  and  Condenser 

(0.0  °C/kW) 

Elevation  cf  Boiler 

15  ft 

Abcwe  Ccncenser 

(4.57  m) 

Channel  Allowance 

0.83  ft 

(0.25  m) 

Efficiencies  for  the  Turbine 

85* 

and  the  Pumps 

Seawater  Temperatures  - Hot 

75°F 

(23 . 9°C) 

- Cold 

4 0°  F 

(4 . 4 °C ) 

Net  Output  Power 

8 ,532x 1C7  BTU/hr 
(25  MH) 

Eundle  Eactor  (F) 

3.0 

Ccst  Indexes  (I) 

1.00  fer  all  ccap. 

Tube  hall  Thickness 

0.028  in 

(0.071  ci) 

Mir.  Distance  Between  Tubes 

0.5  in 

(1.27  cm) 

Heat  Transfer  Coefficients 

145  BTC/f t*-hr-°F 

fer  Eciler  and  Condenser 

(0.822  kii/m2-°C) 

£§t€£ials 

Eciler  and  Ccndenser 

Carbon  Steel  Shells 

Brass  Tubes 

Feed  Pumps 

Bronze 

Circulation  Pumps 

Stainless  Steel 
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Table  6 


EXAMPLE  PEOBLSa 

SOLUTION 

5UHHAFY 

Ja£i3ii€ 

flange 

Mean 

Var 

V 

Units 

Bciler: 

Velocity 

3. 323-3.347 

3.33  5 

0.7 

ft/sec 

(1.013-1.020) 

(1.016) 

(a/sec) 

F ressure 

afFiox  0 

118.2 

0.  C 

lb/in2 

(815.  0) 

(*?) 

lute  length 

39.  16-40  .00 

39.71 

2.1 

ft 

(11.91-12.19) 

(12.10) 

(») 

Shell  cia. 

24.12-25.00 

24.91 

3.5 

ft 

(7.  352-7.620) 

(7.59  3) 

(■) 

lute  dia. 

0. 756-0.761 

0.758 

0.7 

in 

( 1.  919-  1.934) 

(1.926) 

(ca) 

• of  shells 

7.  141-7.409 

7.228 

3.7 

- 

Cct  denser : 

Velocity 

2. 793-2.82 3 

2.808 

1 . 1 

ft/sec 

(0.851-0.860) 

(0.856) 

(»/sec) 

E ressur e 

91.53-91.60 

91.56 

0.1 

lb/in2 

(63  1.  1-631.6) 

(631.3) 

(kp) 

lute  length 

39.35-40  .00 

39.78 

1.6 

ft 

(11.99-12.19) 

(12.12) 

(m) 

Shell  dia. 

44.69-45,00 

44.89 

0.7 

ft 

(13.62-13.72) 

(13.68) 

(») 

lube  dia. 

0.7632-0.7686 

0.7660 

0.7 

in 

(1.938-1.952) 

(1.946) 

(cm) 

1 of  shells 

3. 297-3.401 

3.331 

3.  1 

- 

Hass  flc«  rate 

22.12-22.21 

22.17 

0.4 

10  *lt«/hr 

(10.03-10.07) 

(10.06) 

(10*kg/hr) 

Systea  cost 

7.576-7.673 

7.607 

1.2 

10*  S 
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The  arithmetic  mean  of  th<?  seven  solutions  is  calculated. 
The  difference  between  the  maximum  and  the  minimum  values  is 
shewn  as  a percent  of  the  mean.  The  percent  variaticn  gives 
scie  idea  cf  how  well  the  solution  was  determined.  Note 
that  the  tute  length  and  shell  diameter  constraints  are 
binding  for  both  the  boiler  and  the  ccndenser.  The 
sutprcblea  still  terminates  only  when  the  penalty  function 
does  not  decrease  from  one  iteration  to  the  next.  The 
magritude  cf  ^P(X,r)  is  not  becoming  small  enough  to 
satisfy  the  convergence  criterion.  The  value  never  becomes 
less  thar  C.C2  and  is  usually  greater  than  0.5  when  a 
subpreblea  is  stopped  by  the  program. 

Useful  results  are  obtainable,  however.  The  encineer 
can  learr  much  about  his  problem  even  from  a single  point. 
Fcr  example,  from  the  cost  breakdown  in  Figure  6,  the 
portion  cf  the  costs  attributable  to  each  component  can  be 
seen.  The  heat  exchangers,  of  course,  dominate  with  18%  of 
the  cost.  The  circulation  pumps  are  about  21%  of  the  cost 
with  the  feed  pumps  being  1.2%  of  the  total  cost.  Since 
the  efficiency  assumed  for  the  circulation  pumps  is  prcbatly 
toe  high,  these  pumps  would  be  even  mere  ccstly  in  terms  of 
energy  and  money.  This  could  indicate  the  priorities  that 
should  he  set  on  the  design  efforts  for  improvements  cf  the 
vaiicus  ccmpcnents  with  the  highest  priority  beinc  the 
highest  percentage  cost  component. 

The  system's  power  distribution  is  shown  in  Table  7. 
The  significant  features  are  the  large  energy  absorption 
rate  for  the  siall  useful  output  power  and  the  parasitic 
pumping  pewer.  The  lower  the  efficiency  the  greater  must  be 
the  power  absorbed  for  the  same  output.  This  is  to  be 
expected  fcr  this  kind  of  system.  The  parasitic  pumping 
pewer  (121  cf  the  turbine  output)  is  about  twice  what  it  is 
for  a fcssile  fuel  plant. 


66 


I 


$94,100  PUMPS 

1.2%  $654,000 

8.6% 


FIGURE  8.  COST  BREAKDOWN 


Table  7 

SOLUTION  FOWEF  EUCGET 


Feat  Absorption 
Ttrbice  Shaft  power 
Eeed  Pumps 

Eli  Circulation  Pumps 
Ccnd  Circulation  Pumps 
Feat  Eejected 
Net  tseful  Power 


ETU/hr 

M w a 1 1 s 

3 . 56x  1 C9 

1042 . e 

9 . 6 9x  10  7 

28.4 

4 . 4 9x  1 0 6 

1 . 3 

2 . 58x 1 C6 

0.8 

4 . 5 1 x 1 C6 

1 . 3 

3 . 4 7x  1 C9 

1015.7 

8 . 5 3x  1 0 7 

25.0 
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Cncc  the  mcdel  has  been  optimized  with  the  number  of 
shells  as  real  variables,  the  problem  can  te  re-solved  twice 
by  fixing  tie  shells  at  the  next  higher  and  the  next  lower 
integer  ruaber.  For  the  higher  number  of  shells,  the  cost 
increased  slightly  as  expected,  from  7.607  tillion  dcllars 
average  to  7.983  Billion  dcllars.  The  problem  with  the 
lcwer  integer  number  of  shells  could  not  be  solved  due  tc 
the  constraints  cn  the  tube  length  and  the  diameter  cf  the 
shells  being  binding  before  the  shell  numbers  were 
integerized.  The  program  is  unable  to  find  a new  solution 
that  satisfies  the  eguality  constraints. 

Table  8 shows  some  miscellaneous  data  from  the  mcdel 
that  deserves  comment.  The  overall  heat  transfer 
coefficient  cf  the  bciler  is  approximately  three  tixes  that 
for  the  condenser  which  is  the  major  reason  for  the  lcwer 
ccst  of  the  boiler  relative  to  the  condenser.  However, 
remember  that  the  bundle  factor,  egual  to  3,  is  a part  of 
the  heat  transfer  coefficient.  More  data  must  be  gotten 
abcut  this  factor  in  crder  to  improve  the  reliance  or  any 
sclution.  There  is  a significant  difference  between  the 
Prepane  beat  transfer  coefficients  that  were  assumed  and  the 
theoretical  values  calculated  frem  the  correlations.  Other 
values  should  be  tried  until  the  fixed  value  is  withir  the 
2CJ  tolerance  on  the  theoretical  value.  The  flew  of 
seawater  tbreugh  the  tubes  in  both  the  bciler  anc  the 
ccrdensei  is  turbulent.  Apparently,  the  resulting  cost 
savings  due  to  the  improved  heat  transfer  rate  through  a 
smaller  area  is  more  than  the  increased  costs  of  the  higher 
pumping  pewer  needed.  It  seems  likely  that  the  solution 
obtained  is  globally  optimal;  several  different  starting 
points  were  tried.  The  problem  has  been  started  in  both  the 
laainar  and  the  turbulent  flow  regions  with  no  change  in  the 
results. 
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Table  8 

ADDITIONAL  DATA  FROM  SAHELr 
PROBLEM  SOLUTION 


Boiler  Condenser  Units 
iicrking  Fluid  Heat  Transfer  Coefficients 


Assumed 

145. 

145. 

BTII/h  r-f  1 2-°F 

(0. 822) 

(0.822) 

(kW/a2- cC) 

Theoretical 

227. 

241. 

(1.25) 

(1.37) 

Surface  Area  per  Shell 

2.39x10s 

1.27x10® 

ft2 

(2.22x10*) 

(1.  18x10s) 

(m2) 

Seawater  Cutlet  Temp 

72.3 

43.8 

°F 

(22.4) 

(6.55) 

(°C) 

Shell  Saturation  Teiap 

66.7 

49.9 

°F 

(19.3) 

(9.94) 

(°C) 

Seawater  Pass  Flow  Rate 

4.39x10® 

9 . 11x10* 

lbm/hr 

(1.99x10®) 

(4.  13x  1 0 9) 

(kg/hr) 

Mean  Temp  Difference 

6.67 

7.85 

Of 

(3.71) 

(4.36) 

(°C) 

Overall  Heat  Transfer 

310.6 

105.0 

ETU/hr-f t2-op 

Coef f iciert 

(1.76) 

(0.595) 

(kW/m  2-°C) 

Heat  Flus 

2070.0 

824.9 

ETU/f t2-hr 

(6.528) 

(2.602) 

(kW/m2) 

Seawater  Eeyrclds  to. 

19,800. 

9,960. 

- 
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I V . SECGMHENDATICNS 


The  present  model  has  a great  deal  mere  inforiaticn  to 
offer,  but  scie  improvements  in  the  model  should  be  mace  to 
improve  the  ability  of  a ' nonlinear  program  to  solve  the 
problem.  The  first  recommendation  would  be  to  study  the 
equality  constraints  and  place  inequality  constraints  on  any 
variable  or  grouping  of  variables  if  they  teth  can  tecome 
negative  anc  also  must  have  their  logarithms  computed. 
Examples,  already  included  in  the  program,  are  the  number  of 
tutes  in  tie  shells.  The  heat  transfer  coefficients  for  the 
working  fluid  side  in  the  boiler  and  the  condenser  cculd  be 
included  as  equality  constraints.  This  can  be  dene  by 
adding  the  heat  transfer  coefficients  to  the  vector  I and 
requiring  that  the  variables  be  equal  to  the  coefficierts  as 
computed  ly  the  appropriate  correlation.  The  same  cculd  be 
done  with  Falen  and  Small's  property-dependent  calculation 
of  the  separation  area  by  equation  (12)  and  (13).  From  an 
overall  standpoint,  the  results  seem  to  indicate  either  the 
prcfclem  is  tcc  difficult  to  be  solved  efficiently  ir  its 
present  foci  or  that  another  nonlinear  prcgiam,  other  than 
SORT , might  he  better  suited  to  the  problem.  Numerical 
dif ferentiaticn  of  the  constraints  in  this  problem  do  not 
seem  to  help  the  program  to  find  the  solution  since  the 
magnitude  of  ^ P (X,r)  is  net  near  zero  when  a sukprotlem 
steps . 

Once  the  optimization  scheme  is  functioning  properly, 
the  model  can  te  used  to  study  the  effects  cf  fouling  and 
ccrrosicc  and  the  costs  of  their  prevention.  Many  proposals 
have  keen  made  that  must  be  analyzed  in  the  context  of  the 
total  systei.  Seme  such  as  Amertap,  where  plastic  fcam 
kails  are  pumped  through  the  tubes  to  clean  the  heat 
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transfer  surfaces,  add  a cost  penalty  frci  the  capital  and 
cpe rati  115  ccst  of  the  Amertap  system  and  frcm  the  increased 
pumping  pcwer  required.  Coatings  add  the  capital  ccsts  of 
the  coating  and  the  added  capital  cost  of  increased  heat 
transfer  area  reguired  because  of  the  increased  thermal 
resistance  cf  the  ccating. 

A sensitivity  analysis  can  be  performed  cn  the  model  to 
test  the  assumptions  and  the  various  constants.  The 
sensitivity  tc  changes  in  these  factors  can  indicate  how 
much  one  shculd  be  willing  to  pay  for  the  improvement 
sought.  high  sensitivity  indicates  areas  where  the  greatest 
effort  fcr  improvement  should  te  desired  or  the  priorities 
shculd  te  set  cn  research  ictc  the  unknown  areas  cf  the 
problem.  The  trade-off  between  cost  and  performance  can  be 
studied  by  changing  the  output  power  and  looking  fcr  the 
lowest  ccst  per  kilowatt-hour,  which  is  the  measure  by  which 
ccrventicral  pcwer  plants  are  compared.  One  could  lcok  into 
the  cost  and  the  benefits  of  increasing  the  depth  frcm  which 
the  cold  water  is  pumped  versus  the  improved  temperature 
differential  possible  from  usicg  the  colder  water. 

For  the  capital  cost  of  the  system  to  be  complete,  the 
turbine  design  must  be  included  and  ccst  estimates  made. 
The  turbine  will  be  unique  due  to  the  very  low  pressure 
potential  drop  it  must  use.  Some  work  has  been  dore  by 
ether  greups  working  on  OT£C,  Hef.  25  and  26.  Reference  27, 
by  Bcbert  1.  Bartlett,  can  provide  costing  and  design 
information  although  he  deals  only  with  steam  turbines.  The 
remaining  part  of  the  system  that  is  desireable  for 
inclusicc  ir  the  analysis  presented  herein  is  the  cold 
seawater  intake  pipe.  It  will  be  on  the  order  of  1,0C0  ft 
(3C4.6a)  lerg  and  100  ft  (30.4  m)  in  diameter. 
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